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Abstract 

By means of the quark mass density- and temperature- dependent model, it 
is found that the negative charge and the higher strangeness fraction are in 
favor of the stability of strange quark matter at finite temperature. A critical 
baryon number density hbc has been found. The strange quark matter is 
unstable when < ns c - The charge and strangeness fraction dependence of 
stability for strange quark matter is addressed. 
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I. INTRODUCTION 



It has been found from the theoretical study that the introduction of strange quarks into 
a plasma with two flavors could lower the Fermi energy of the system and thus the mass 
of quark matter Jl] . If its mass is lower than the mass of hyperonic matter with the same 
strangeness fraction, or even lower than the non-strange nucleonic matter, the strange quark 
matter (SQM) would be the true ground state of QCD. Therefore, to search the SQM in the 
laboratory is an important task for recent heavy ion collision experiments at ultrarelativistic 
energy @. 

The essential problem for detectability of SQM in heavy ion experiments is to study its 
stability during the formation of quark gluon plasma (QGP). Many properties of stability 
of SQM, for example, the "stability window" [H-H], the charge-dependence and the 

strangeness fraction- dependence p3|-|T4|j of stability have been widely discussed by many 
authors. They employed different models to investigate this problem. But unfortunately, 
many conclusions are different, in particular, part of them, say, the charge dependence is 
contrary. 

The electric charge is of vital importance to the experimental searches of SQM. By 
means of MIT bag model, Jaffe and his co-workers ||||,|7j suggested that the SQM is slightly 
positively electric charged. They considered strong and weak decay by nucleon and hyperon 
emission together and concluded that the stable SQM will have a low but positive charge to 
mass ratio. Using the same model as that of Jaffe et al. but considering the initial condition 
of possible strangelet production in relativistic heavy ion collision carefully, Greiner and his 
co-workers [0,|l2| argued that strangelets are most likely highly negatively charged. They also 
emphasized that the stability of infinite SQM depends on the bag constant B. An infinite 
SQM, treated as a gas of noninteracting quarks, is absolutely stable for B l l A ~ 145MeV, and 
metastable for B 1 ^ ~ 150 — 200MeV, but unstable for larger bag constants. Of course, this 
conflict of conclusions is limited in the framework of MIT bag model and at zero temperature. 

Employing another model, namely, the quark mass density-dependent (QMDD) model, 
Peng and his co-workers reconsidered this problem |J. They introduced a critical density 
and proved that SQM cannot maintain its flavor equilibrium, and thus no longer be stable 
under this critical density. They found that proper negative charges can lower the critical 
density and thus make the SQM become more stable. Their conclusion is not surprised if 
one notice the confinement mechanism of QMDD model. In QMDD model, the mass of u, d 
and s quarks are given by [|15H T7|,3,|9|,[7|J^1 

m q = - — , (q = u,d,u,d), (1) 
3n B 

B 

m S! s = m s0 + - — , (2) 
3n B 

where n B is the baryon number density, m s o is the current mass of the strange quark matter 
and B is the vacuum energy density inside bag (bag constant). The mechanism of con- 
finement can be mimicked through the hypothesis Eqs.(p]) and (Q). The masses of quarks 
become infinitely large as the volume increases to infinity or the density decreases to zero. 
So that the vacuum is unable to support it j|,[17). Noting that the boundary condition 
of confinement of MIT bag model corresponds to the zero quark mass inside the bag but 
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infinity at the boundary or outside the bag, we see that the confinement mechanisms are 
similar to these two models. Due to this correspondence, as pointed out by Benrenuto and 
Lugones Q], in almost all cases, the properties of SQM given by QMDD model are nearly 
the same as those obtained in the MIT bag model. The conclusion of Peng et al. supports 
the result given by Greiner et al. 

This paper evolves from an attempt to extend these discussions to finite temperature. 
The reasons of this extension are: at first, the SQM might be possible and formed during the 
phase transition from hadronic matter to a deconfined QGP. From lattice QCD calculations, 
the phase transition of deconfinement will happen at critical temperature. Secondly, as 
was pointed out by our previous papers PJIBfl, QMDD model cannot be used to explain 
the process of the quark deconfinement phase transition because the quark confinement is 
permanent in this model. Instead of this model, a quark mass density- and temperature- 
dependent (QMDTD) model in which the quark confinement is not permanent has been 
suggested by us p|,|TB|j. As its application, we hope to employ QMDTD model to discuss the 
charge and strangeness fraction dependences of SQM at finite temperature. 

The organization of this paper is as follows: In next section, we review the QMDTD 
model and give the main formulae which are necessary for studying the thermodynamical 
behaviors of SQM. In section 3, we give the result of numerical calculations for QMDTD 
model and investigate the stability of hot SQM in detail. Our conclusions are summarized 
in the last section. 



II. QUARK MASS DENSITY- AND TEMPERATURE-DEPENDENT MODELS 



To overcome the difficulties of QMDD model which appear at finite temperature, in refs. 
PJT8|], based on the Friedberg-Lee model we suggested a QMDTD model. The basic 
difference between our model and QMDD model is that, instead of a constant B in Eqs.(|l]) 
and @, we argued that B must be a function of temperature and choose B(T) as 



B = B Q 
B = 0, T > T c , 



'-0 



, < T < T c 



(3) 
(4) 



where Bq is bag constant at zero temperature and T c is the critical temperature of quarks 
deconfinement. Substituting Eq.@ into Eqs.(|l]) and (|), we find the masses of u, d, s quarks 
and their corresponding anti-quarks as 



Bp 
3n B 



-'0 



(q = u, d, u, d), 



m s0 + 



3n B 



rp x 2 



(5) 
(6) 



when < T < T c . When T > T c , the masses of quarks are independent of temperature: 
m u,d,u,<i — and m St s = m s o- Due to the difference between Eqs.(|l]),(@) and Eqs.©,(|]), as 
was proved in ref. 0, we found that when n B approaches to zero, the temperature T — > T c 
in our model, but diverges in QMDD model. In QMDTD model, the masses of u, d quarks 
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become zero (or constant for strange quarks) when the temperature approaches T c . This 
means that the quarks can be deconfined. Here, the vacuum energy density B is a decreasing 
function as temperature increases. As was pointed out by Greiner et al., the chosen values 
of B affect the stability of SQM remarkably ||. Therefore, it is of interest to consider the 
effect of B(T). 

Treating the infinite SQM as a gas of noninteracting quarks and electrons HH, we find 
that the thermodynamical potential density of SQM is 

g(T r°° dNi 

(2^ Jo dk 

where i stands for u, d, s (or u,d,s) and the electron e(e + ), gi = 6 for quarks and antiquarks, 



" = E «< = - E £L f m (i + .-*.»-«>) , (7) 



gi = 2 for e and e + . £i(k) = ym 2 + k 2 is the single particle energy and rrii is given by 
Eqs.(|^),(^j) for quarks and antiquarks. /i; is the chemical potential (for antiparticle \i\ = 
—fj>i). The number density of each particle can be obtained from eq.(|7j) by means of 



i on 

V dpi 



(8) 



and we find 



n- r°° / 1 1 



Am = m — m = - — — / d k — ^ tttt n • (9) 

For SQM, the baryon number density satisfies 

n B = ^(An u + An d + An s ). (10) 

The electric charge density Q of SQM reads 

2 11 
Q = -An u - -An d - -An s - An e , (11) 

and the light quarks are converted to strange quarks through the weak processes P,|16| 

u + d < — > u + s, s — > u + e~ + u e , d — > u + e~ + i> e , u + e~ — > d + u e . (12) 



Neglecting the contribution of neutrinos ^,|], W , 17 , we obtain the conditions of chemical 
equilibrium from the reactions ([12]) as 

A* s = Vd, (13) 

Hs= Vu + He,- (14) 

Noting that the chemical potentials of particles and anti-particles are determined by 
Fermi distribution Eq.(|), we see that Eqs.©, (|), (|), (0), (0), ® and (0) form an 
equation group and can be solved self-consistently. The solution of this equation group 
determines the chemical stability of SQM. The results of numerical calculation will show in 
the next section. Here we want to point out that at a fixed temperature T and for a fixed 
value of charge density Q this equation group which determines the configuration of SQM 
has no solution under a critical baryon number density ns c - This means that the SQM 
is unstable when Ub < n.B C . With the values of Q changed and the temperatures limited 
between < T < T c , the condition n# > ub c determines a stable area which might be useful 
for the detectability of SQM in experiments of high energy heavy ion collisions. 
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III. CONCLUSION AND DISCUSSION 



The numerical calculations of the equation group mentioned above have been done with 
the parameters set: 

B = 170MeVfnT 3 , m s0 = 150MeV, T c = 170MeV. (15) 

Firstly, we study a SQM system with neutrality charge Q = 0, for a given n B , the 
fractions of different quarks are given by 

Fq(q=u,d,s) = (16) 

At a fixed temperature T = lOOMeV, the u, d, s quarks fraction are drawn in Fig.l by solid 
line, dashed line and dotted line, respectively. At high baryon number density, all of the 
u, d, s quarks tend to become a triplicate. When the density becomes lower, u fraction 
increases and s fraction decreases monotonously. We see from this figure that the fractions 
of u, d, s quarks satisfy an inequality: 

F d >F u > F s . (17) 

In particular, we would like to emphasize that, as shown in Fig.l, three different quark 
fraction curves stop at a critical baryon number density (n Bc = 8.28 x 10~ 3 fm~ 3 ). One can 
not find real roots for F d , F u and F s provided n B < n Bc . At n Bc , F s = 0.107 ^ 0. In ref. ||, 
the authors defined the critical density by the condition F s = 0. Of course, the SQM can 
not maintain the chemical equilibrium if no strange quarks exist. But as shown in Fig.l, 
this is only a sufficient but is not the necessary condition. Because, one cannot find solution 
in the ranges < F s < 0.107. 

To compare the QMDD model in which B is a constant and our QMDTD model in which 
B satisfies Eqs.(||) and (f|), we plot the critical density n Bc vs temperature T with neutral 
charge density Q = in figure 2, where the dashed line and the solid line refer to QMDD 
model and QMDTD model, respectively. We see that the curve given by QMDD model can 
never be zero even the temperature increases to T > T c . It can easily be understood if we 
notice that m q approaches to infinity when n B —>■ from Eqs.(|I]) and (0). To excite a particle 
with infinite mass must spend infinite energy and then infinite temperature. Because of 
n Bc 7^ 0, when n B < n Bc , the system would still be unstable even at very high temperature 
for QMDD model. But for QMDTD model, we see from the solid line that n Bc = at 
T = T c = 170MeV. It means that the SQM is completely stable and may be detected from 
high energy heavy ion collision. 

Now we are in a position to address the charge dependence of stability for SQM. The 
critical density vs temperature curves with various charges are shown in figure 3, where the 
solid line, dotted line and dashed line refer to Q = 0, — O.lefm -3 , — 0.2efm~ 3 , respectively. 
We see that the critical density n Bc decreases rapidly and the stable area expands widely 
when Q decreases to more negative. So the negative charge is in favor of the stability of 
SQM. Our result supports the conclusion given by Greiner et al. || and extends their result 
to finite temperature. 

To study the effect of the strangeness fraction, we show the same curves as that of Fig.l 
but for Q = — 0.2efm -3 in figure 4. Comparing Fig.l with Fig. 4, we find that all of three 
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curves tend to become a triplicate for Q = — 0.2efm -3 at high density, too. And there 
are two major differences between these two figures, firstly, the critical density decreases 
from 8.28 x 10~ 3 fm~ 3 to 5.02 x f -4 fm~ 3 when charge density Q changes from neutrality to 
negativity — 0.2efm~ 3 ; secondly, Eq.(|TTD becomes 

F d > F s > F u . (18) 

The fraction of s quarks is higher than the fraction of u quarks. This is of course reasonable 
because s quark has negative charge and u quark positive. The critical density n Bc decreases 
as the F s increases. The higher strangeness fraction is in favor of the stability of SQM. To 
illustrate this result more transparently, using the expression of strangeness fraction 

/, = ^, (19) 
n B 

which is usually employed by many references, we draw the critical strangeness fraction f sc 
vs temperature T curves in figure 5, where the solid line and dashed line refer to Q = and 
— 0.2efm -3 , respectively. The f sc is strangeness fraction corresponding to the critical density 
riBc at the same conditions of charge and temperature. We see from Fig. 5 the critical 
strangeness fraction increases with increasing temperature. The stable area corresponds to 
the larger negative charge and higher strangeness fraction. 

Finally, we hope to investigate the stability of SQM from another view point. Instead 
of considering the chemical equilibrium, we study the behavior of chemical potentials for 
electron and u quark, respectively. The reason for this choice is that the electron (u quark) 
is the particle with minimum mass but negative (positive) charge in our model. The smaller 
the particle mass is, the lower the chemical potential will be. When chemical potential tends 
to be negative, the system becomes unstable because the Fermi surface is negative in this 
case. Therefore, we can use the conditions \x e > 0, \i u > to determine the stability of 
SQM. 

The chemical potential of electron \i e vs baryon number density %ata fixed temperature 
T = lOOMeV with various positive charge density Q — 1 x 1CT 4 , 1 x 1CT 3 , 3 x 10~ 3 , 5 x 1CT 3 
and 7 x 10~ 3 efm -3 are drawn in figure 6. We see that fj, e becomes negative and the system 
becomes unstable when Q < O.OOlefm -3 . Similarly, the chemical potential of u quark /i u 
vs baryon number density at a fixed temperature T = lOOMeV with various negative 
charge density Q = —0.2, —0.322, —0.4, —0.6 and —0.8 efm -3 are drawn in figure 7. We see 
from this figure that when Q < —0.322, the curve will pass through a negative area and 
thus SQM is unstable. Therefore, at temperature T = lOOMeV, the charge density must be 
limited to O.OOlefm" 3 > Q > -0.322efm" 3 for stable SQM. 



IV. SUMMARY 



In summary, by means of the QMDTD model, we find that the negative charge and the 
higher strangeness fraction are in favor of the stability of SQM at finite temperature. This 
result is in agreement with that given by Greiner et al. || at zero temperature, but we extend 
their conclusion to finite temperature. A critical baryon number density n^ c under which the 
SQM is unstable has been obtained. We find that ns c decreases rapidly when charge density 
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Q becomes more negative and approaches to zero at the deconfinement temperature T c . By 
using the chemical potentials of electron and u quark, we suggest a method to determine 
the limit of charge density for the stability of SQM. We hope the stable areas given by our 
model can have impact on the detectability of SQM in future heavy ion experiments. 
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VI. FIGURE CAPTIONS 

Fig.l The quark fraction F q as a function of baryon number density ub with neutral 
charge, F u , F d and F s are represented by solid, dashed and dotted lines, respectively. Here 
F d >F u > F s . 

Fig. 2 The critical density n B as a function of temperature T with neutral charge, the 
dashed line is for QMDD model and the solid line is for QMDTD model. 

Fig. 3 The critical density % as a function of temperature T with various charge den- 
sity Q = 0, —0.1 and — 0.2efm~ 3 , they are represented by solid, dotted and dashed lines, 
respectively. 

Fig. 4 The quark fraction F q as a function of baryon number density ub with a negative 
charge density Q = — 0.2efm~ 3 , F u , F d and F s are represented by solid, dashed and dotted 
lines, respectively. Now F d > F s > F u . 

Fig. 5 The critical strangeness fraction f s as a function of temperature T. The solid line 
is for Q = and the dashed line is for Q = — 0.2efm~ 3 . 
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Fig. 6 The chemical potential of electron function of baryon number density n B 

with the positive charge density Q — 1 x 1CT 4 , 1 x 1CT 3 , 3 x 1CT 3 , 5 x 10~ 3 and 7 x 10~ 3 efuT 3 . 
At a fixed temperature T = lOOMeV, they are represented by solid lines respectively. 

Fig. 7 The chemical potential of u quark fj, u as a function of baryon number density n B 
with the negative charge density Q = —0.2, —0.322, —0.4, —0.6, and — 0.8efm~ 3 . At a fixed 
temperature T = lOOMeV, they are represented by solid lines respectively. 
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